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Abstract 

Half-lives ^'^'^ neutrino positron double beta decay modes /3'^EC/ECEC are calculated 
for ®^Sr, a nucleus of current experimental interest, within the framework of the deformed shell 
model based on Hartree-Fock states employing a modified Kuo interaction in (^^3/2, ^/5/2i ^Pi/2j 
^59/2) space. For a reasonable description of the spectra of *^Sr and ^^Kr and to generate allowed 
GT strengths, the single particle energies of the proton and neutron ^^9/2 orbitals, relative to the 
^P3/2 orbital energy, are chosen to be (3.5 MeV, 1.5 MeV) for both '^^Sr and ®^Rb and (1.5 MeV, 1.5 
MeV) for ^^Kr. With this, the calculated half-lives for the /3+EC and ECEC modes are ~ lO^^yr 
and ~ 4 X lO'^^yr respectively. 

PACS numbers: 23.40.Hc, 21.10.Tg, 21.60.Jz, 27.50. +e 
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I. INTRODUCTION 



Double beta decay (DBD) is a rare weak interaction process in which two identical nucle- 
ons inside the nucleus undergo decay either with the emission of two neutrinos or without 
any neutrinos. The two neutrino double beta decay {2vl3~ which was first predicted 
long back by Meyer [1] has been observed experimentally in more than 10 nuclei [2, 3]. In 
contrast, the positron decay modes, i.e. 2u /3+/3+//3+EC/ECEC decay modes (hereafter, 
all these three combined is called 2v e"'"DBD) are not yet observed experimentally (excep- 
tion being ^^°Ba decay derived from geochemical methods [4]). However, in the last few 
years there are serious attempts, using direct counting methods, to measure half-lives for 
2v e+DBD modes in many nuclei ranging from ^^Zn to ^^^Ba [5-9]. In the mass A=80 
region, the candidate nuclei for 2v e"'"DBD are ^^Kr, ^^Se and ^^Sr. For ^^Kr and ^^Se there 
are already some experimental efforts and following this, recently [10, 11] we have carried 
out calculations for these two nuclei using the so called deformed shell model (DSM). In 
addition we had previously attempted to study 2i//3~/3~ nuclear matrix elements for ^^Ge 

^^Se using DSM [12]. There was also an attempt for ^^Se ^^Kr and it is seen that [13] 
DSM produces nuclear matrix elements within a factor of 2 the QRPA results given in [14]. 

Over the years, the DSM model based on Hartree-Fock states has been used to study with 
considerable success: (i) spectroscopic properties, such as band structures, shapes, nature of 
band crossings, electromagnetic transition probabilities and so on in A=64-80 nuclei [15-20]; 
(ii) T — 1 and T — bands in N=Z odd-odd nuclei and T — 1/2 bands in odd-A nuclei by 
including isospin projection [21-23]; (iii) transition matrix elements for // — e conversion in 
^^Ge [24] and in the analysis of data for inelastic scattering of electrons from /p-shell nuclei 
[25]; (iv) /3-decay half hves, GT distributions and electron capture rates in N ~ Z nuclei 
in A=60-80 region [10] and so on. It is seen that the predictions of DSM for 2i> e"'"DBD 
half-lives for ''^Kr are close to those of QRPA and PHFB models. For ''"^Se there are no 
other model predictions except those of DSM. See [10, 11] for details. 

A group led by H.J. Kim at Korean KPNU facility has initiated DBD experiments for 
positron double beta decay (both 2u and Ou) in ^^Sr using SrCl2 crystals [26, 27]. Prompted 
by these efforts, we have carried out DSM calculations, extending the studies in [10, 11] 
further, for 2z/ e"^DBD half-lives for ®'^Sr nucleus and the results are reported in this paper. 
It should be pointed out that there are no theoretical results available in literature for this 
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nucleus. Just as with "^Kr and " Se, we have undertaken the 2v e+DBD study of ^ Sr as 
it was aheady shown that DSM model describes very well [28] the low-lying bands in this 
nucleus. Also as we shall discuss ahead, the proton and neutron single particle energies (spe) 
of the ^(?9/2 orbital play an important role in describing the spectroscopic results of ^'^Kr and 
in determining the ®^Sr ^^Kr 2u e"'"DBD half lives. Now we will give a preview. 

In Section II given is, for completeness, a brief discussion of DSM model and the formalism 
for calculating 2v e+DBD half-hves. Spectroscopic results for ^^Sr and ^^Kr are given in 
Section III. In Section IV, DSM results for 2v e+DBD half lives are discussed and also given 
here are some concluding remarks. 



II. 2v e+DBD HALF-LIVES AND DEFORMED SHELL MODEL 

Half-life for the 2i/ e+DBD decay modes for the Oj 0^ transitions, is given by [29], 

[T!^k)Y' = G^.{k) \M2.f (1) 

where k denotes the modes P'^P'^, and ECEC. The kinematical factors G2v{k) are 

independent of nuclear structure and they can be calculated with good accuracy [29, 30]. 
On the other hand, the nuclear transition matrix elements (NTME) are nuclear model 
dependent and they are given by [29, 31], 

where |0+), |0+) and |l+> are the initial, final and virtual intermediate states respec- 
tively and En{Ej) is the energy of the intermediate (initial) nucleus. Similarly, Eq — 
|Wo where Wq is the total energy released for different 2i/ e''"DBD modes. As given 
in [31-33], Wo{/3+/3+) = g^+^+ + 2me, Wo{/3+EC) = Q^+eq + and lyo(ECEC) = 
^ECEC ~ 2?7T.e + e^i + ■ The Q- values for different 21/ e+DBD modes arc [31]: 
Q^+^+ = M{A, Z) - M(A, Z -2) - 4me, Q^+ec = M{A, Z) - M{A, Z - 2) - 2m^ and 
Qecec = M{A, Z) — M{A, Z — 2) . Here M denotes the neutral atomic mass (available for 
example from the tabulations in [34]) and is the binding energy of the captured atomic 
electron. Energies in the denominator in Eq. (2) are taken in the units of electron mass. 
With the atomic mass difference being 1787 keV for ^^Sr decay, it should be clear that the 
/3'^/3'^ DBD mode is forbidden (as the Q-value will be negative) for ^^Sr. For ^+EC and 



ECEC modes, one and two K shell electron capture respectively are considered. The binding 
energy for the K shell electrons is taken from [35] . 

In DSM, for a given nucleus, starting with a model space consisting of a given set of 
single particle (sp) orbitals and effective two-body Hamiltonian, the lowest prolate and 
oblate intrinsic states are obtained by solving the Hartree-Fock (HF) single particle equa- 
tion self-consistently Excited intrinsic configurations are obtained by making particle-hole 
excitations over the lowest intrinsic state. These intrinsic states will not have good angular 
momentum and hence good angular momentum states are obtained by angular momentum 
projection from these intrinsic states. The normalized states of good angular momentum 
projected from the intrinsic state Xk{v) can be written in the form 

^mAv) = ^^^^ J dQD-I^Amm \xk{v)) (3) 

where Njk is the normalization constant given by 

2 14-1 1"^ 

Njk = d/3sm/3diM {xK{v)\e-'^'^\xK{v)) ■ (4) 

In Eq. (3) fl represents the Euler angles {a, (5, 7), R{Q) which is equal to 
exp{—iaJz)ex.p{—if3Jy)ex.p{ —ijJz) represents the general rotation operator. In general the 
projected states with same J but coming from different intrinsic states will not be orthog- 
onal to each other. Hence they are orthonormalized and then band mixing calculations are 
performed. DSM is well estabhshed to be a successful model for transitional nuclei (with 
A=64-80) when sufficiently large number of intrinsic states are included in the band mixing 
calculations. Performing DSM calculations for the parent, daughter and the intermediate 
odd-odd nucleus (here we need only the 1+ states) and then using the DSM wavef unctions, 
the (TT~ matrix elements in Eq. (2) are calculated. For further details of DSM see [10]. 

Let us add that the recently introduced projected configuration interaction (PCI) model of 
Horoi et al [36, 37] and the projected shell model (PSM) of Sun et al [38-40] are quite similar 
to DSM (PCI may be better for quasi-spherical nuclei and PSM includes BCS correlations 
from the beginning). 
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III. SPECTROSCOPIC RESULTS FOR ***SR AND ^^''KR 
A. ^Sr 

In the DSM calculations for the structure of the nuclei ^"^Sr, ^^Rb and ^^Kr and for 
^^Sr 2u e'^DBD half-lives presented in section IV, we have used a modified Kuo effective 
interaction [41] in the (^^3/2, ^/5/2, ^fi'9/2) space with ^^Ni as the inert core. The 

spe of these orbitals are taken as 0.0, 0.78, 1.08 and 3.5 MeV respectively in a first set 
of calculations and in the next set of calculations, the ^5*9/2 energy for the protons and/or 
neutrons is changed to 1.5 MeV for reasons discussed in detail ahead. DSM with modified 
Kuo effective interaction has been quite successfully used by us in describing many important 
features of nuclei in A ~ 60-80 region. In particular, shape coexistence in spectra, observed 
B{E2) values, identical bands, band crossings, and deformations in 76,78,80,82,84gj. jgotopes 
are well described by DSM [17, 28]. First we will describe briefly the earlier spectroscopic 
results for ^^Sr. 

Starting with the modified Kuo interaction and ^^3/2, ^/5/2, ^Pi/2, ^99/2 spe (both for 
protons and neutrons) to be 0.0, 0.78, 1.08 and 3.5 MeV respectively, DSM calculations 
have been carried out for ^^Sr in [28]. Here, first axially symmetric HF calculations are 
performed and the lowest prolate HF intrinsic state is obtained. The reason for neglecting 
the oblate states has been discussed in detail in [19]. The lowest HF sp spectrum for ^^Sr 
is shown in Fig. la; note that the HF sp states with a given k values will be doubly 
degenerate. By particle-hole excitations from the lowest intrinsic state shown in Fig. la, 
excited intrinsic states are generated. In the band mixing calculations a limited number 
of five configurations are taken and they are: (i) ground — 0'^ intrinsic state; (ii) two 
excited K'^ — 0"^ intrinsic states obtained by exciting two or four protons into the k states 
arising out of the ^99/2 orbit; (iii) a = 8'^ intrinsic state obtained by exciting two valence 
neutrons into k = 7/2^ and 9/2+ states; (iv) a K'^ = 2"*" intrinsic state obtained by exciting 
a valence proton from k = 5/2" into 1/2^ state. As discussed in detail in [28], the band 
mixing calculations reproduce the observed data quite well. Figure 2 shows the calculated 
and experimental spectra. The DSM results in Fig. 2 are shghtly different from those in Fig. 
2a of [28] where the so called tagged HF calculation is done (in the present calculations, just 
as in [10, 11], no tagging has been done). The results for 10"*" and 12+ levels are in better 
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agreement with data compared to the earher results in [28]. The two lowest excited 8"^ bands 
in Fig. 2 are estabhshed by experiments [42, 43] to be two neutron and two proton ahgned 
bands respectively and DSM reproduces this structure. Similarly, the observed B(E2) values 
and also the excited 2+ band are well described. All these confirm that DSM gives good 
spectroscopic results for ^"'Sr. 

In the 2v e+DBD decay, the intermediate nucleus is ^^Rb. Fig. lb gives the HF spectrum 
for ^^Rb (for completeness we also show the HF spectrum for ^^Kr in Fig. Ic). As we can 
see from Fig. 1, in the DBD process, in the lowest order, the valence proton from the 5/2~ 
orbit in ^^Sr will go to the unoccupied neutron 7/2+ orbit in ^^Rb. This gives a negative 
parity state for ^"^Rb and this is not allowed by GT. Hence an intrinsic state with a proton 
in 5/2+ in *^Sr is needed to generate allowed GT matrix elements. However the intensity 
of this configuration will be small in the ®^Sr ground state. Thus the GT matrix elements 
will be small for ^^Sr to ^^Rb. Therefore we need to reduce the neutron spe for these 
two nuclei so that the negative parity neutron sp states will be near the Fermi surface. As a 
result we will be generating many low-lying K'^ — 1+ levels in ^^Rb and there will be large 
GT matrix elements for ^^Sr to ^"^Rb (similarly also for ^'^Rb to ^^Kr). Following this, we 
have considered proton ^5(9/2 spe e'^^^{}g%/2) = 3.5 MeV and neutron ^g()/2 spe e*-"^("'^5'9/2) = 1-5 
MeV. The HF sp spectrum for this choice is shown in Fig. 3a. Starting with the lowest 
configuration shown in Fig. 3a, five more excited configurations are considered in the band 
mixing calculations as before for ^^Sr. Fig. 2 shows that practically there is no effect of 
the change in £^"'^(^^'9/2) on the spectrum (also wavefunctions and hence the results in [28] 
are well preserved) . The reason being that there is essentially no change in the proton and 
neutron occupancies. In Fig. 2, deviations of the results of the earlier calculation from the 
present calculation are shown in the parentheses with an asterisk. Deviations below 50 keV 
are not shown. 

In conclusion, for ^^Sr and ^^Rb nuclei, for DSM calculations we employ the spe for ^^3/2, 
^/5/2 and ^pi/2 orbitals (same spe for both protons and neutrons) to be 0.0, 0.78, 1.08 MeV 
respectively while [6*^^^(^(79/2), £("^(^5(9/2)] = [3.5, 1.5] MeV. Now we will consider ^^Kr where 
there were no earlier DSM results. 
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B. s^Kr 



Employing the same set of spe as chosen for ^^Sr and ^ Rb, the lowest prolate HF intrinsic 
state is obtained for ^^Kr; note that as above e^^'^ (^519/2) = 3.5 MeV and e^") (^519/2) = 
1.5 MeV. Then the HF sp spectrum in Fig. Ic remains unchanged except that the HF 
neutron 15^9/2 fc-states move down in energy by 2 MeV. By particle- hole excitations from 
the lowest intrinsic state, excited intrinsic states arc generated for ^^Kr and band mixing 
calculations are carried out using lowest 15 configurations. The observed lowest 6"*", 8"*" 
and 10+ states are two-proton aligned states and similarly the 8^ and 10^ are two-neutron 
aligned bands. It is seen that the calculated two-proton aligned band is high (by more than 
1 MeV) compared to experiment. Also the 6+ level produces weak B{E2) for 6^*" — > 4]*'; the 
calculated value is 0.2 W.u. while 6.9 W.u. is the data value (see Table 1). In order to bring 
the two-proton aligned band head to be close to data, it is necessary to lower the e^^\^g9/2) 
energy to 1.5 MeV in addition to using ^^^^^99/2) = 1-5 MeV. The HF sp spectrum with 
[e^)(^g'9/2), e*-"^(^5'9/2)] — [1-5, 1.5] MeV is shown in Fig. 3c and the energy spectrum is shown 
in Fig. 4. In the DSM calculations, starting with the lowest configuration shown in Fig. 3c, 
15 excited configurations are considered for band mixing. It is seen from Fig. 4 that the 
two-proton aligned band is in reasonable agreement with experiment. From Fig. 4 and also 
the B(E2) results in Table 1, wc conclude that DSM gives a reasonably good description 
of the spectroscopy of ^"^Kr with the choice [e*-*'-'(^5'9/2), e*-"-'("'^5'9/2)] = [1-5, 1.5] MeV. Let us 
add that the recent data on occupancies in ^^Ge and ^^Se do point out that [44, 45] lg'9/2 
is probably much closer to the fp shell than anticipated in the past, thus justifying the 
lowering of the energies of proton and neutron I519/2 orbitals. It is also pluasible that part 
of the reason for these lower energies may be due to the neglect of excitations into l5'7/2 and 
2(i5/2 orbitals. 

In conclusion, for ^^Kr for DSM calculations we employ the spe for V3/2, ^/5/2 and ^pi/2 
orbitals (same spe for both protons and neutrons) to be 0.0, 0.78, 1.08 MeV respectively 
while [e(f)(V9/2),e(")(V9/2)] = [1-5,1.5] MeV . 
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IV. RESULTS FOR 2v e+DBD HALF LIVES AND CONCLUSIONS 



Using the DSM wavef unctions, 2v ^+EC and ECEC half-lives for ^^Se ^^Kr transitions 
are calculated and the results are shown in Table 2. In the DSM calculations for the inter- 
mediate nucleus ^^Rb, 24 — 1^ bands are mixed and they span 3 MeV from the lowest 
K'" = 1+ band. The resulting 24 1"*" states are employed in Eq. (2) to calculate the GT 
matrix elements. Wc have verified that the inclusion of additional l"*" states do not change 
significantly the half-life results. Further, the 0"*" ground state of ^"^Sr is generated by mixing 
30 intrinsic states. Similarly for generating the 0"^ ground state of ^^Kr 11 intrinsic states 
are used in the DSM calculations. The Kinematical factors are taken from [30]. It is 
seen from Table 2 that the predicted half- fives are ~ lO^^yr and this should be of interest 
for future experiments. Instead of the choice [e^'^^{}g^/2)-,e^'^''{}g9/2)\ = [1-5, 1.5] MeV, if we 
use (3.5, 1.5)MeV for generating the ground state 0^ of ^^Kr just as the spe used for (^^Sr, 
^■^Rb), the half-lives reduce by a factor ~ 2 as shown in the brackets in Table 2. Let us add 
that the results in Table 2 are the first nuclear structure results for 2v /3+EC and ECEC 
half-fives for «^Se ^ ^^Kr. 

In conclusion, prompted by recent experimental interest in ^^Sr 2v e+DBD and the results 
reported in [10, 11] for '''^Kr and ^^Se using the DSM model, we have carried out spectroscopic 
calculations and then the DBD half-lives calculations for ^^Sr — )■ ^^Kr. We have shown that 
the DSM model gives reasonably good spectroscopy for both ^''Sr and *^Kr by adjusting the 
spe of proton and neutron orbitals. Proceeding further we have presented the DSM 
results for half-lives for ^+EC and ECEC modes. These results should be a good guide for 
the experiments initiated by Kim et al [26, 27] at KPNU, Korea. Finally, formulation within 
DSM, for calculating Oi/DBD and e+DBD half lives is being developed and the results 
for ^"^Sr e+DBD half lives will be reported in future. 
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TABLE L B(E2) values for ^^Kr in Weisskopf units (W.u.). DSM results are compared with 
experimental results (called EXPT in the table) from [43]. The effective charges employed are 
Cp = 1.6e and e„ = l.Oe just as in the previous paper on ^^Sr [28]. 



Transition DSM EXPT 



2^- 




18.4 


12 


4^- 


.2+ 


28.9 


15 


6|- 


.4+ 


18.8 


6.9 


2^- 


^2+ 


3.2 


13.3 


<- 


^2+ 


4 


1.6 


n- 


.0+ 


1.6 


3 



TABLE II: Deformed Shell Model results for half-lives T^J^ and the corresponding phase space 
factors for possible decay modes for ®*^Sr — ^^^Kr. Note that the natural abundance of ^^Sr 
is 0.56% [46]. The G21/ and half- lives are calculated using QA/gv = 1- For half-lives given in the 

brackets see text. 

Mode G2u (in yr"^) T^J^ (in yr) from DSM 

/3+EC 1.504 X 10-2^ 1.2 x 10^^ (5.41 x lO^^) 
ECEC 4.367 x 10-^3 4.16 x 10^^ (1.87 x lO^^) 
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FIG. 1: HF single particle spectra for ^^Sr, ^^Rb and ^^Kr for e(P)(^59/2) = e^"H^59/2) = 3.5 MeV. 
In the figures circles represent protons and crosses represent neutrons. The Hartree-Fock energy 
(E) in MeV, mass quadrupole moment (Q) in units of the square of the oscillator length parameter 
and the total K quantum number of the lowest intrinsic states are given in the figure. 
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FIG. 2: (a) Calculated and (b) experimental spectra for ^^Sr. Experimental data are from [42, 
43]. Note that the DSM spectrum is obtained using [e^P\^gg/2),e^'^\^g9/2)] = [3.5, 1.5]MeV. The 
numbers in the brackets (over a few levels) indicate the difference between the calculated energies 
and those obtained with the choice 99/2) 1^^"^^^ 99/2)] = [3.5, 3.5]MeV. See text for details. 
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FIG. 3: HF single particle spectra for ^^Sr, ^^Rb and ^^Kr. Note that [e'-P\^gg/2),e^''\^gg/2)] = 
[3.5, 1.5]MeV for ^^Sr and ^^Rb nuclei and [e^^'^C 99/2), (^^""H^ 99/2)] = [1-5, 1.5]MeV for ^^Kr. In the 
figures circles represent protons and crosses represent neutrons. The Hartree-Fock energy (E) in 
MeV, mass quadrupole moment {Q) in units of the square of the oscillator length parameter and 
the total K quantum number of the lowest intrinsic states are given in the figure. See text for 
details. 
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FIG. 4: (a) Calculated and (b) experimental spectra for ^^Kr. The experimental data are from 
[43]. See text for details. 
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